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Abstract 

By defining dijet rapidity gap events according to interjet energy flow, we 
treat the photoproduction cross section of two high transverse momentum jets 
with a large intermediate rapidity region as a factorizable quantity in perturba- 
tive QCD. We show that logarithms of soft gluon energy in the interjet region 
can be resummed to all orders in perturbation theory. The resummed cross 
section depends on the eigenvalues of a set of soft anomalous dimension matri- 
ces, specific to each underlying partonic process, and on the decomposition of 
the scattering according to the possible patterns of hard color flow. We present 
a detailed discussion of both. Finally, we evaluate numerically the gap cross 
section and gap fraction and compare the results with ZEUS data. In the limit 
of low gap energy, good agreement with experiment is obtained. 

PACS Nos.: 12.38.Aw, 12.38.Cy, 13.85.-t, 13.87.-a 
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1 Introduction 



In a recent paper we have presented an explanation, from perturbative QCD, of 
the dijet rapidity gap events experimentally observed in pp scattering at the Fermilab 
Tevatron 0. These events consist of a pair of jets produced at very large momentum 
transfer, and separated by a wide empty region, where barely any particle multiplicity 
is measured. They were originally predicted from the exchange of at least two gluons 
in a color singlet configuration, which avoids color recombination between the jets 
and associated interjet hadronization 0. A quantitative formalism to describe gap 
events, however, has been lacking in the literature. 

In Ref. [0 we have discussed the dependence of the dijet cross section on the energy 
flow into the central region, clearly related to the observed particle multiplicity. We 
have shown how to factorize the partonic cross section into a hard scattering function, 
accounting for the dynamics of the highly virtual quanta exchanged in the scattering, 
and a soft function, describing the soft radiation emitted into the interjet region. 
Both functions are defined as matrices in the space of the possible color flows at 
the hard scattering, formalizing the idea, already expressed by other authors, that 
the color structure of the short-distance partonic scattering is not uniquely defined 
f|, ||. The soft function contains, at each order in perturbation theory, logarithms 
of the interjet radiation energy, which can be resummed, to give the dependence of 
the dijet cross section on the interjet energy flow. The resummation is driven by 
the eigenvalues and eigenvectors of a soft anomalous dimension matrix, defined in 
the space of color flows. In Ref. 0, we have performed the analysis of this matrix 
for qq — ► qq scattering, which is the relevant partonic process in a valence quark 
approximation for pp collisions. The result indicates that, in the limit of a very large 
interjet region (corresponding to a very large parton center of mass energy, compared 
to the momentum transfer), the color singlet component dominates, thus merging 
our picture, in this asymptotic configuration, with approaches to the problem based 
on the Regge limit of QCD |J and on the idea of color singlet dominance || [F]. 
We have already pointed out that a full treatment of the Tevatron dijet gap cross 
section requires us to include the contribution of processes involving gluons, which 
we postpone, for that specific problem, to forthcoming work ||. 

In this paper, we will apply our method to analyze the photoproduction of two 
high transverse momentum jets with a large intermediate rapidity gap. Events of this 
kind have been reported in e + p collisions at HERA |J . Here, the quantity of interest 
is the gap fraction, the ratio of the number of dijet events with a large interjet gap, 
to the total number of dijet events. The ZEUS experiment shows that, binning the 
gap events according to the width of the interjet empty region, the gap fraction, after 
an initial decay, reaches an asymptotic plateau. It has been shown in Refs. |6| that 
the fraction of gap events from the exchange of a color octet, reggeized gluon falls off 
exponentially with increasing gap width. The left-over asymptotic excess has been 
interpreted as the fraction of the color singlet exchange of a hard BFKL pomeron |5J. 
We emphasize here that the BFKL approach is complementary to ours, because it 
deals with the resummation of gluonic radiation, which, in our language, would be 
described as hard (see Sec. ^ below). 
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Applying the methods of Ref. |]J to this problem is especially interesting, because, 
while the data for the Tevatron gaps show the dependence of the cross section on the 
interjet energy at fixed gap width, here we will explore the opposite dependence, 
variable gap width at fixed energy We will show that the factorization and resum- 
mation of the soft gluonic energy emitted into the gap leads precisely to the behavior 
of the gap fraction observed by the experiment. In addition, photoproduction is an 
ideal process in which to analyze the contribution of gluons to gap cross sections. 
Theoretical studies of jet photoproduction have been pursued since the early 80's 



10 , p/Lj |12| , |T3| , [14], [15[ and are a very active area of research. The dominant partonic 



mechanisms in photoproduction are the direct and resolved scatterings: in the former 
case the photon from the incoming lepton interacts directly with the quark or gluon 
from the proton, and in the latter it fluctuates first into a hadronic state of low vir- 
tually, acting as a source of partons, which then scatter off partons from the proton. 
The precise experimental determination of the partonic content of the photon is still 
an open problem, partly because it requires isolating the resolved component of the 
scattering |L(], |j"6f . However, existing data indicate that, especially at low values of 
partonic x, the dominant component of the photon is gluonic fL2| . Correspondingly, 
we will see that for the resolved contribution, the scattering of a gluon from the pho- 
ton with a quark from the proton is dominant in the kinematical region of interest. 
Also the contribution from the resolved reaction g + g —>■ g + g has to be taken into 
account, although it turns out to be much less important in this case, because of the 
high transverse energy cuts experimentally imposed on the jets. 

One last important remark has to be made, about the issue of survival R ITTA. In 



dijet gap events, the survival probability is limited by the probability of no radiation 
inside the gap from the interaction of spectator partons. Such probability is estimated 
to be of the order of 10% in pp scattering, if we require a truly empty interjet region 



H \F7j. Presumably, for ep reactions the number of spectator quarks and gluons 
is reduced with respect to pp scattering, because of the partly leptonic nature of 
the initial state. Therefore our resummed formula, which accounts for the amount 
of perturbative survival from soft gluon dynamics, should be less sensitive to these 
effects. 

Throughout the paper, we will work in close correspondence with the ZEUS exper- 
imental configuration described in Ref. f9j. In Sec. |], we will discuss the kinematics of 
the problem, define the dijet gap cross section, and review its factorization properties. 
In Sec. D, we will introduce the gap fraction, as the ratio of the dijet gap cross section 
to the overall dijet cross section. We will also identify the partonic reactions giving 
the largest contributions to both the gap and the overall cross section. For each of 
these reactions, in Sec. (|, we will present the explicit decomposition into hard and 
soft parts. In Sec. [5], we will present the process-dependent soft anomalous dimension 
matrices, their eigenvalues and eigenvectors, which govern the soft dynamics of the 
scattering. Finally, in Sec. we will give numerical results for the overall dijet cross 
section, the dijet gap cross section and the gap fraction. We will draw a comparison 
with the experimental results of Ref. M, and present our conclusions. 



3 



2 Dijet Gap Cross Sections in Photoproduction 



2.1 Definitions 

In this section we will introduce dijet gap cross sections in photoproduction. We first 
recall the definition of particle pseudorapidity, 



V 



where 9 is the angle of a particle momentum with respect to a fixed direction, typically 
the beam direction. We will consider positron-proton scattering 




e + {p A )+p{p B ) — - e + (p' A ) + Ji(pi) + J 2 {P2) + X, 



Sap- 



(2) 



for the production of two jets at fixed pseudorapidity difference, Arj — r\\ — r\i. We 
sum inclusively over final states, while measuring the (soft) energy flow into the 
intermediate region between the jets. The boost-invariant pseudorapidity difference 
A?7 fixes the partonic scattering angle 9 through the formula 



— — = In (cot ( — 
2 I I 2. 



(3) 



Following Ref. |J, the jets are defined by cones of radius R = 1.0 in the rj-cp plane 
mapping the ZEUS calorimeter detector. This constrains the interjet region to have 
(pseudo)rapidity width 

Ay = At? - 2R. (4) 

The low-virtuality photons exchanged at the electromagnetic vertex of the inter- 
action can be thought of as real particles of energy E n = yE e , whose spectrum is 
given by the Weiszacker- Williams formula |18L |i~9| 



F j/ e (y) 



l + (l-y) 2 ln (QL x (l-y) 



2tt 



y 



mly 2 



+ 2m e y 



(i-y) 



Q 



2 2^ i2 

max '"'eel 



(5) 



Here a em is the electromagnetic coupling, m e is the electron mass, and Q max * s the 
maximum photon virtuality, determined from the (anti) tagging conditions of the 
experiment. From Ref. we take the value Q max = 4GeV 2 . The generic differential 
cross section for electron-proton scattering can then be viewed as a convolution of 
the photon distribution in the electron and the photon-proton cross section: 



da ep (S) 



dyF y/e (y) da yp (S^ p ), 



(6) 



where S is the center of mass energy squared of the electron-proton system, and where 
5* 7 p = yS is its analog for the photon-proton system. We use for y min and y max the 



values given in Ref. M, y a 



0.2 and y n 



0.8. 
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2.2 The factorized cross section 



For our specific case, the inclusive cross section for dijet events with transverse energy 
greater than Et, rapidity difference At], and energy flow less than Q c in the interme- 
diate region, of rapidity width Ay, can be written as follows [[10], [□], [TJ, [L3|, [14]]: 

CPr " c A ^ " r dyF l/e {y) 



dAr] 



(Q c ,S,E T ,Ay) 



2/min 

dAr] 



(Q 



Et, Ay) + — 2. (Q c 



dAr] 



S 1P , E T , Ay) 



(7) 



Here da^/dAr] and da^/dAr], are the direct and resolved contributions to the cross 
section respectively. In the former case the low-virtuality photon interacts directly 
with the parton from the proton, while in the latter, as mentioned above, it acts itself 
as a source of quarks and gluons, which then scatter off the partons from the proton. 
The two cross sections can be written in factorized form as [|K], 11], [12, [13], [14] 



da d[T 



S 1P , Et, Ay) 



E 

fp,fl,fr 

X dAr] 



dx p (f>f p / p (x p , -i) 



Qc,i,s,r]jj,Ay,a s (i) 



and 
da 



dAr] 



S<yp, Eq 



Ay) 



E 



fix 



dx p 0/ 7 / 7 (x 7 , -t) <pf p /p{x p , -i) 



x [Qc t, s, rjjj, Ay, a s {t) 



(9) 



In these formulas 0/ 7 / 7 and <pf p / p are parton distributions in the photon and proton 
respectively, evaluated at the scale —t, the dijet momentum transfer, which is related 
to the partonic center of mass energy squared, s, and the dijet rapidity difference, 
Ar], according to the formula i = — | (l — tanh(^)). The integration regions for the 

partonic fractions x p and x 7 are denoted by R Xp and R Xl respectively. da^/dAr] 
and do^/dAr] are hard scattering functions, starting from the lowest order Born 
cross section. The index f (7f) denotes the partonic process / 7 + f p — > f\ + fi 
(l+fp fi+h)- The detector geometry constrains the phase space for the dijet total 
pseudorapidity, r]jj = {j]j 1 +r]j 2 )/2, with \r]jj\ < 0.75 ||. Similarly, the lower bound 
on the transverse energy of the jets, Et, and the dijet pseudorapidity determine the 
phase space for the partonic center of mass energy squared s, with 4E T 2 cosh 2 
s < S yp exp(2r]jj)yE e /E p . 



< 



2.3 Refactorization of the partonic scattering 

In this section we will review some of the arguments already presented in Refs. 
|H|, 20], [H, 22| , to show how it is possible to perform a further factorization on the 
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partonic scattering functions, da^'/dArj and da^/dArj, of Eqs. (§) and (|9|). The un- 
derlying argument is that in the partonic scattering the soft gluon emission decouples 
from the dynamics of the hard scattering, and can be approximated by an effective 
cross section, in which each of the partons is treated as a recoilless source of gluonic 
radiation. Formally, this is equivalent to replacing each parton with a path-ordered 
exponential of the gluon field (eikonal or Wilson line) in the proper representation 
of SU(3), the fundamental representation 3 (3*) for quarks (antiquarks), and the 
adjoint representation for gluons. In this way the hard amplitude is replaced by a 
sum of eikonal operators, u>/(x){ Cfc }, depending on color tensors cj, which account 
for the color flow at the hard scattering, times short- distance coefficient functions. 
The effective dimensionless eikonal cross section ( in the following f will refer to both 
direct and resolved processes) can be written as 



x(o|f [K(o)){ M ] IO(^Mo) {6i} ] |o> 



(10) 



where we sum over all the final states subject to the constraint of having energy 
less than Q c in the interjet region of rapidity Ay, while cutting off all the remaining 
integrals at the ultraviolet scale /i. This makes the eikonal cross section free from 
potential collinear singularities associated with gluon emission from the Wilson lines 
H, |21|, |22j. The Latin indexes / and L refer to the color structures of the amplitude 
and of its complex conjugate. At tree level, with no soft gluons, the above formula 
reduces to the square of eikonal vertices, with matrix elements given by traces of the 
color tensors in the amplitudes. 

In these terms, the partonic scattering function can now be factorized into the 
product, in the space of color tensors, of a hard scattering matrix [] , H®, accounting 
for the quanta of high virtuality exchanged in the scattering, and a soft matrix, S L j, 

da® , , 

-j^j (Qc,t,s,rijj,Ay,a s (t)) = —1^2 cosh 



xHfl (AvW^ 2 )) fiS Uv,-) , (ii) 




where we follow the convention of the sum over repeated indices. This factorization 
holds to leading power in -k-, with A the QCD scale parameter. We can identify a 



hard scale, y — t, a soft scale, Q c , and a new factorization scale, /i. The soft matrix, 
Sli, precisely coincides with the gauge invariant eikonal cross section of Eq. (|T0"1), 




a& eik) (-^) • (I 2 ) 




1 The normalization of the dimensionless hard scattering matrix, Hj^, can be found in Eq. (|68|) 



IL 

below. It differs slightly from the one used in Ref. 
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In general, corrections to the factorized expression in Eq. (JTTj) are expected from 
three-jet final states, but, since the cost of adding an extra jet to the final state is at 
least one power of a s (t), they will be suppressed. 



2.4 Evolution of the soft matrix in color space 

The left-hand side of the factorized expression, Eq. (0), is independent of \i. This 



means that the /i-dependence of the two matrices, Hfl and S^j, must cancel in the 
product. We thus derive for S^j the evolution equation [|] 



7 (f) 
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dfx 



s 



(f) 

LI 



fp( f M c( f ) ow (r^>\ 

I 1 S ILB^BI ~ ^LAK 1 S )AI 



(f) / r ( f )> 



(13) 



where (a s ) is a process-dependent soft anomalous dimension matrix. Solving this 
equation will enable us to resum all the leading logarithms of the soft scale Q c . It 
is convenient to treat Eq. (|T3| ) in a basis which diagonalizes r^(a s ). Following Ref. 



|22|| , we denote by Greek indexes the basis of the eigenvectors of T 



«.), {i4 f) )}> 

corresponding to the eigenvalues Xg . We then transform Eq. ( |TTD to this basis and 
solve the evolution equation (|i~3"f) for by integrating with respect to /x between 

the soft scale Q c and the hard scale y—i, to get: 



ci<r (f) 

[Q c ,s, i, rjjj, Ay, a s (i) 



71 



— ( 2 cosh' 
2s 



v 2 



x E H f« (Ay, v 7 !, V-f, a s (t) ^ 0) (Ay) 



hi 



Qc 
A 



(f) 




(f) 



(14) 

The double differential cross section of Ref. 0, d 2 a^ /dAr]dQ c , giving the distribution 
of dijet events as a function of interjet radiation and rapidity interval, is obtained by 
differentiation with respect to Q c . In the above formula the kinematical cuts on the 

jets require the minimum value of to be exactly Et (Et = 5GeV in Ref. ||), 
still much larger than A. The exponents E a L are given by 



E<£1 (Ar/, Ay) 



2tt 



X^*(Ar,,Ay)+Xf (Ar], Ay) 



(15) 



where (3\ is the first coefficient in the expansion of the QCD /3-function, [3\ = ^N c — 



]rif, and where we define XP by A 



a s X ( P + 



Observe that we have explicitly 



mantained a distinction between the dependence on Ar] and Ay, in spite of the fact 
that for this specific problem they are related through the identity in Eq. (f|). In 
fact At? can be seen as a parameter of the hard scattering (related to the scattering 
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angle), whereas Ay, at fixed Arj, is a measure of the width of the interjet region and, 
as such, has a geometrical interpretation. The matrix S^k in Eq. (|14]) is obtained by 

transforming to the new basis the zeroth order S^f^ of Eq. (|TTJ) , which, as mentioned 
above, is just a set of color traces. The change of basis goes through the matrix 
((i?( f )) -1 )^ = ( e /j J^j according to the formula p2 



c-(f.o) 



((^ (f) )- 1 )1 ^S((J2W)-%. (16) 



r(f.i)> 



Analogously, we take for the if}^ s the squares of the Born- level amplitudes, rep- 
resented in the original color basis. In the diagonal basis, each hard matrix H\l 
becomes Ha^ , defined by 

= (17) 

Observe that S {i ' 0) and if (fil) both acquire a Ay- dependence through the change of 
basis. 

The relevant A7/-dependence in the right hand side of Eq. (O) is contained in the 
exponents E a p, Eq. dl5|). We will see below that, for all the partonic processes, the 
real parts of the eigenvalues of Vp (a s ) in the exponents are positive definite functions, 
with boundary value Re fA^) = at Arj = 2 (Ay = 0). Then, since Q c < y—t, 
the last line in Eq. (lj) acts as a suppression factor. Specifically, the eigenvalue with 
the biggest real part will drive the strongest suppression for the component of the 
scattering oriented along the direction of the corresponding eigenvector. As we shall 
see, most of the relevant information about the At] behavior of the cross section is 
derived from the properties of tl 
we will study in detail in Sec. §. 



derived from the properties of the soft anomalous dimension matrix rV(ai s ), which 



3 The Gap Fraction 

We will be interested in the evaluation of the gap fraction, defined as the ratio of the 
number of dijet events with a specified rapidity gap to the total number of dijet events. 
A gap event is usually identified experimentally by the lack of particle multiplicity 
in the interjet region [Q, [)]. The multiplicity is determined from the number of 
calorimeter cells which measure energy deposition above a threshold, 300 MeV in the 
case of Ref. |J. Therefore, the condition for a gap event is the absence of such cells 
(or clusters of cells) in the rapidity region between the jets. Our formulation of the 
problem, Eq. (|7|), is in terms of the total interjet flow, Q c , of hadronic radiation, which 
is clearly related to particle multiplicity. In analogy with experiment, we introduce an 
energy threshold, Qq, which is different in principle from the experimental calorimeter 
threshold, and identify a gap event from the condition of interjet radiation less than 
Q . From Eqs. (|7]), @ we see that the dependence on the interjet radiation is all 
in the partonic scattering, Eq. flTj|). Then, in order to get a dijet gap cross section, 
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3.0 



Figure 1: The overall dijet cross section compared with the experimental data of Ref. 



we just need to evaluate Eq. (|TJ|) at the threshold value Q of the interjet energy flow 

Qci 

i 



^(f),gap 

^ (Qa, s, i, T]jj, Ay, a s (t) 



71 

2§ 



2 cosh' 



X 



In 



Qo 
A 



?( { ) r 



In 



A 



The denominator of the fraction is also obtained from Eq. ([H]), this time evaluated 

at the maximum value of Q c , the hard scale \f—i. Values of Q c above this scale 
would break the premise underlying the factorization in Eq. (|TT|), since the emission 
into the interjet region would not be soft any more. By performing the substitution 

Qc = m Eq. (|14D we find the full leading order partonic dijet cross section 



(M, ^j, Ay, a s (t) 



71 



2 cosh' 



Hl 1] Uy, Vs, V^, a s (ij) S^ 0) (Ay) , 



(19) 
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2 

Ar| 



Figure 2: The contribution of the different partonic reactions to the dijet cross section. 
At Arj = we can identify from top to bottom: the overall result (solid line); the 
contributions of: 7 + g — > q+q (dotted line) , q + g — > q + g (dashed line), q + g — > g + q 
(dot dashed line), g + g — > g + g (double dot dashed line), 7 + q(q) — > g + g(g) (short 
dashed line), g + q — > g + g (short dotted line), g + g — > g + g (short dot dashed line). 
Here the two reactions q + g — > q + g and g + g — * g + g differ from each other by 
the exchange of the Mandelstam invariants t and u (see Eqs. (|69|) and (|70|) in the 
Appendix). 



where, from Eqs. (0) and (JH]), for each partonic process the trace identity i?£^, S^f 1 

Hj f £^ S^IP holds, which is proportional to the leading order partonic cross section (see 
also Sec. |] and Appendix |A|). After summing the gap and total cross sections over 
partonic subprocesses, as in Eqs. (|8|) and (^), the gap fraction has the form 

^(Qo,S,E T ,Ay) 

/■gap _ dAr] v^ u ' ' 1 ' »> 

d ^(S,E T ,Ay) 

In the numerical evaluation of the denominator of this fraction, we use, for each 
partonic reaction f in Eq. (|19D , the invariant S-matrix elements squared, and corre- 
sponding parton luminosities, summarized for completeness in Appendix |A|. We 
have used four quark flavors, rif = 4, and have assumed flavor symmetry for the sea 
quarks. The numerical result we have found, shown in Fig. [T], is in reasonable agree- 
ment with the data points from Ref. [Q. Here it should be emphasized that, for the 
overall dijet cross section, and also, as we will see in Sec. || for the gap cross section, 
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corrections are expected from next-to-leading order contributions to the hard matrix 



flOl , ITT] , [L2| |T^| . However, the gap fraction should not be too sensitive to them. 

In Fig. |2] we show the partonic reactions most relevant to the overall cross section, 
evaluated on the full Ar^-range. For the A^-range of Ref. |J, A77 > 2.0, Fig. |3] 
clearly shows that it is sufficient to consider the contributions from the direct process, 
7 + g — > q + q, the resolved gluon-initiated reactions, q + g — > q + g and g + g — > g + g, 
and the quark processes, q + q — > q + q and q + q — > q + q. These will also be 
the dominant contributions to the numerator of the fraction, which differs from the 
denominator only by the last two factors in Eq. (|18D , depending on the soft dynamics 
of each specific partonic reaction. The determination of the gap cross section will 
require an analysis of the partonic scattering in terms of color flow, to be found in 
the next Section, and the knowledge of the spectrum of the process-dependent soft 
anomalous dimension matrices, to be discussed in Sec. ||. 



4 Lowest Order Hard and Soft Matrices in Color 
Space 

In this section, we present the color decomposition of the hard scattering for the 
resolved partonic processes that have been shown above to give the dominant con- 
tribution to the dijet cross section in the At] region of interest. For each case, we 
will define a basis of color tensors, ||c^/)|, and give, in this basis, the Born level 

hard matrix, H\l , whose elements come from the squares of the color-decomposed 
tree amplitudes. We will also give the zeroth order soft matrix, S^[f \ defined by the 



set of traces S^'P = Tr 



The knowledge of these matrices is not only 

crucial for the evaluation of Eqs. ( |TBD and (|I13D, but will be important for threshold 
resummation in jet production as well 



For the color tensors, following Ref. [£4] , we will find it convenient to use t-channel 



bases. In terms of momentum and color, we label the partonic process f according 
to / 7 (l A , r A ) + f p (l B , r B ) -> fi (pi,n) + f 2 O2, r 2 ), where {l A} r A }, {h, r B } ({pi,n}, 
{P2, ^2}) are, respectively, the momenta and colors of the incoming (outgoing) partons. 
The hard matrix will be expressed in terms of the partonic Mandelstam invariants 

s = (l A + l B f 
i = (l A -pi) 2 

u = (l A -p 2 ) 2 ■ (21) 

It will also depend on the QCD running coupling evaluated at the hard scale, (— 1) 1//2 , 
according to the formula 

/~\ 27T 

ot,[t) = 7 ; r- 22 

W ftln^-^Va/A) 

For each process, it is straightforward to check that the color trace over the product 
of hard and soft matrices gives, as expected, the tree level invariant matrix element 
squared of Appendix |A] |24 . 



11 



3.0 



Figure 3: The dijet cross section (solid line) and the result obtained by dropping from 
Fig. the contributions of the direct process 7 + q(q) — > g + q(q), and of the resolved 
reaction q + g — > g + q (dotted line). 



4.1 Hard and soft matrix for qg — ► qg 

For the partonic process q + g—>q + gwe define the basis of color tensors 



d 

ifr 



(23) 



Cl 

c 2 

c 3 = ifr B r 2 c(Tp) rirA , 

where Ci is the t-channel singlet tensor, c<i and C3 the symmetric and antisymmetric 
octet respectively. Here and below, we suppress color indices on the c/s. In this 
basis, the Born level hard scattering can be described in color space by the matrix 

/ 5X1 §Xi fX2 \ 
# (D V§, a,(t)) = — a 2 (t) \ X i Xi 2 X 2 

V 5X2 2 X2 X3 / 

with Xi, X2 and X3 functions of the partonic Mandelstam invariants: 

i 2 

2 — 

su 



(24) 



Xi 



, S M 2 1 t 2 
X2 = 1 - J - - - — 

£ st 2 su 
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su i 2 . , 

Xs = 6-8-^---. 25 
t 2 su 

It is easy to see that, in the limit of forward scattering (At] — > oo), the component of 
the hard matrix in the color direction of the antisymmetric octet becomes dominant, 
since it is 0(i~ 2 ). 

The zeroth order soft matrix, on the other hand, is given by 



S (0) 



( N C (N 2 - 1) 

^- (iV 2-4)(iV c 2 -l) 1. (26) 

\N C (N 2 - 1) 



where N c = 3 is the number of colors. 

4.2 Hard and soft matrix for gg — > gg 

For this process, a suitable t— channel basis of color tensors has been defined in Ref. 



23 to be: 



1 

C\ = ~ [fr A r B ld ri r2l d rA r B lf rir2 l] , 
1 

c 2 = T [fr A r B ld ri r 2 l d r A r B l j ri^i] j 

^3 J [fr A rild rB r2l d rAri lf rBr2 j\ , 

c 4 = Pi(r A , r B ; n, r 2 ) = ^ ri4ri 5 rjjr2 , 



C5 = P8 s ( r A,rB;r 1 ,r 2 ) = -d rAriC d rBr2CJ 



c 6 — P8 A (. r A, r B] r l, r 2) — 7;fr A ricfr B r2C 



3 
— ( 

5 
1 

3" 



C 7 — -PlO+To( r ^' r B'i f\i r 2) = 2 (^ A r s ^rir 2 — ^r A r 2 ^r B ri) 

^^fr A rxcfr B r2C , 

C8 = Prj{ r Ai r B] r l, r 2) = 2 (^r B ^ir 2 + ^r A r 2 ^T- B n ) 

1 3 

g < ^ r -4 ri ^ r s r '2 ~^d rAriC d rBr2C . (27) 

The last five elements of the basis are the t-channel SU(3) projectors for the decom- 
position into irreducible representations of the direct product 8® 8, which corresponds 
to the color content of a set of two gluons. 

In this basis, the hard matrix has the block-diagonal structure 
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where the 5x5 submatrix is given by 



H, 



(i) 

5x5 



-i,Vs~,a s {t^J = Yg a s(£) 



/ 


9xi 


|xi 


|X2 





-3xi 


\ 




Q 

ixi 


Ixi 


|X2 





-fxi 








|X2 


X3 





-1X2 























V 


-3xi 


-|xi 


-fX2 





Xi 


) 



(29) 



with xii X2 and %3 defined by 



tu st 



t 1 



Xi 

X2 
X3 



1-TT7-T77 + — 



St 
~2 



U 



iu u 2 
s 2 st 



su 
tu 



27 (su ltu 1st 



2 V tii 2 sm , 



(30) 



Once again, it can be noticed that the component of the hard scattering in the color 
direction corresponding to the antisymmetric octet dominates in the forward limit. 

For the zeroth order soft matrix, straightforward color traces with N c = 3 colors 
give 



S ( ) 



( 


-5 




















\ 







-5 




























-5 




























1 




























8 




























8 




























20 





\ 























27 / 



(31) 



We notice that the first three color tensors of the basis, Eq. (|27|) , have negative 
eigenvalues, while the eigenvalues of the projectors count the number of color states 
belonging to each irreducible representation. In fact, ci, c 2 and c 3 of Eq. (|27|) de- 
couple from the physical cross section. Eq. (p8|) supports this interpretation, since 



the components of the Born level hard scattering along these color directions vanish 
identically, reducing the dimensionality of the problem, in color space, from 8x8 
to 5 x 5. In principle, from Eqs. (0) and (0), the original dimensionality could be 
restored after the change from the color basis to the basis of the eigenvectors of the 
soft anomalous dimension matrix, I" '5. We will show below that the structure of the 
specific Ts for g + g —>■ g + g prevents this from happening. 



4.3 Hard and soft matrix for qq — > qq and qq —> qq 

Both processes can be conveniently treated in the t— channel singlet-octet basis |22 
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In color space, the hard matrix has the form 



(32 1 
^X2 : 



81 Xl 9X2 
X3 



(33) 



where, in the case of qq — > qq, \\i X2 and %3 are defined by 



Xi 



i 2 + u 2 
~1 2 ~ 

h 2 8i 2 + u 2 



X2 = »Tr-o 



st 3 s 2 
J 2 + u 2 2t 2 + u 2 4« 2 
* 3 = 2 ^ + 9^"3^' (34) 

while in the case of qq — > qq, related to qq — > qq by the crossing transformation s <-> u, 
they are given by 

i 2 + s 2 

Xi = — 1^— 

tr 



X2 



s 2 8s 2 + 1 2 



til 3 u 



2 



s 2 + m 2 2s 2 + t 2 4 s 2 
t 2 9 u z 3tu 

In both cases, for forward scattering, the component of the hard matrix along the 
color octet direction becomes dominant. 

For the lowest order soft matrix, after the color trace, we find the common simple 
form |L| 

5 The Soft Anomalous Dimension Matrix 

In this section we shall describe how to compute the process-dependent soft anomalous 

dimension matrix, 1^ |§ . The //-independence of the left hand side of Eq. (|TT1) implies 

(f) (f) 

that the matrices H}£ and S L j renormalize multiplicatively, 



H ,m = ( Z m-i) 40 [ ( zft r . u 

S<'»S» = (zf% B S^A, (37) 



where the superscript (B) identifies the bare quantities, and Z SCD is a matrix of 
renormalization constants, describing the renormalization of the soft function, i. e. , 



15 




A 1 1 A 



1—2 B 





B 2 — 2 B 




B 2 










Figure 4: Real corrections to the eikonal scattering, Eq. (|i~0l) . For brevity, we only 
show half of the contributing cut diagrams. The remaining ones can be obtained by 
hermitian conjugation of each of the above graphs, i. e. , by reflection with respect 
to the final state cut. 



the eikonal scattering, Eq. flT2p. The one-loop anomalous dimension is obtained 

(i) 

from the residue of the UV pole contained in the matrix Z s , 

( r ") lf (9) = -f| Res -»( 4 ")« (9 - E) ' (38) 

From Eq. ([10|) we see that the potential sources of UV divergences are the virtual 
vertex corrections to the eikonal color- dependent operators wj(x)^ Ck }, and also the 
real corrections, when gluons are emitted into the forward region of the scattering, 
because the eikonal cross section, Eq. (P~D|) , is completely inclusive with respect to the 
amount of forward radiation. 

The Feynman rules for the evaluation of eikonal diagrams have been presented in 
Ref. P2| . The necessary one- loop calculations are illustrated in Figs. given the 
two color vertices cj and cl, corresponding to the hard amplitude and to its complex 
conjugate, we build the eikonal cross section from the product of the color- dependent 

eikonal operators, wi(x){ Ck y and (wi(x){ C| .}j , Eq. flT0), and we compute the diagrams 



obtained by adding one extra real or virtual gluon, shown respectively in Fig. |] and 
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A 1 




B 



Figure 5: Virtual corrections 
show half of the contributing 
hermitian conjugation of eacl 
to the final state cut. 



1 A A 1 




2 B B 2 



1 A A 1 




2 B B 2 



1 A A 1 




2 B B 2 



to the eikonal scattering, Eq. 
cut diagrams. The remainin 
ti of the above graphs, i. e. , 



1 A 




2 B 



(|ToD . For brevity, we only 
g ones can be obtained by 
by reflection with respect 



H]. Working in Feynman gauge, the self-energy diagrams do not contribute, because 
they are proportional to the invariant mass of each lightlike eikonal line. 

From the loop integration, it is easy to see that the UV divergences from the 
gluon emission diagrams and from the virtual diagrams with the gluon emitted and 
reabsorbed between an incoming and an outgoing eikonal line, are real. On the 
other hand, the coefficient of the UV pole from virtual diagrams with the gluon 
emitted and reabsorbed between two incoming or two outgoing eikonal lines, also 
have imaginary parts. Since the imaginary part of the soft anomalous dimension 
matrix comes only from virtual eikonal diagrams, it can be extracted directly from 
Ref. where contributions to T$ are computed, for most partonic processes, from 
the virtual corrections to the eikonal operators wj(x)^ Ck y. The real parts of the UV 
divergences coming, respectively, from virtual and real emission diagrams partially 
cancel in the sum, leaving a remainder which is a function of Ay, the width of the 
interjet rapidity region, where the gluon radiation is measured, and of the scattering- 
dependent parameter Arj. 

We will now present, for each of the partonic processes relevant to the dijet cross 
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(f) 

section, the result for T s , and the analysis of its eigenvalues and eigenvectors, which 



we will relate to the behavior of the cross section, Eq. ([14}). For the direct reaction 
7 + g — > q + q, which has only color octet content, T s will reduce to a function, 



related by crossing to the soft anomalous dimension computed in Ref. p5| for direct 
photon production. 

We will see below, for all the resolved reactions, that the eigenvectors of Ts are 
mixed color states, and that their color composition is a function of the gap width 
only, through the "geometrical" parameter Ay. We have already pointed out in Eq. 
(|j) that, in the experimental configuration of Ref. ||, Arj fixes Ay, leaving us with 
only one independent parameter. We should emphasize, however, that for the case 
of the Tevatron rapidity gaps Jl|, ||, where the gap is taken fixed in the calorimeter 
detector, and is independent of the dynamics of the scattering, the eigenvectors have 
the same form as for this problem. Below, we will find it convenient to refer to the 
eigenvectors as "quasi-color" states (quasi-octet, quasi-singlet, and so on), identified 
by their behavior in the limit of asymptotically large gap regions, where they reduce 
to pure color states. We will discover that the eigenvalue with the smallest real 
part always corresponds to the quasi-singlet eigenvector, and we will show, from the 
discussion in Sec. |2.3| , that the related component of the hard scattering is the only 
one to survive for large values of Ay. 



5.1 Soft anomalous dimension for qg —> qg 

In the basis of Eq. (^), the one- loop soft anomalous dimension matrix is 



T 5 {Art, Ay) 



, Pi + & 

«s 



2tt 



4m 





Pi 



-Z7T 



2m 

N c m 

Pi 



where the functions p\ and £i are defined as follows: 
^(Ay) = N c (i7r-2Ay), 



Pi{Aq,Ay) 
_Nl + l 

2N r . 



2 In 



tanh(^) + tanh(^) 
v tanh(f2) - tanh (^) 

in[ tanh(^) + tanh(^) ^ ^Jl 



v tanh(^) - tanh 



+ m 



2Nl 



tanh (^) 



1 + tanh 



(39) 



(40) 



While the first function depends only on the rapidity width of the interjet region, the 
second also depends on the hard scattering through At]. 
The eigenvalues of the above matrix are given by 



\i{Ar),Ay) 



2m 

UAr,,Ay) = ^ 



h e 3 si + s 2 



2pi + £l / 4i7T 



e 3 s 1 + s 2 
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Figure 6: Plot of the real (left) and imaginary (right) part of the eigenvalues (Aj = 
Xi/a s ) of the soft anomalous dimension matrix for q + g — > q + g scattering. The solid 
line identifies the quasi-singlet eigenvalue, Ai, the dashed and short dashed lines the 
two quasi-octets, A 2 and A 3 . 



A 3 (A77,Ay) 



2tt 



2pi + e 



1 / 4i7T 

- + ( si + e 3 



S-2 



where si and s 2 , in the Ay range of Ref. || 0, Ay < 2, are the functions 



si{Ay) = 



3 7T 2 



^ + ^(^-8)6+ ^(^-4)^+ 



2T4 6 2 (iV c 4 -28iV c 2 + 88)+^(iV c 2 + 4) J 



27 



(41) 



(42) 



and 



s 2 (Ay) = 



27 

2ttJ 
27 



7T 



7T 



27 



28iV c 2 + 88) + !L(iV 6 2 + 4) 3 



(43) 



The eigenvectors corresponding to the eigenvalues of Eq. fl41|) can be written in 



2 The definitions of si and S2 given in Eqs. ( |42| ) and (O) arc valid for Ay < 4tt^^-, the branch 
point of the square root on the right hand side of Eqs. ( f!2[ ) and ( f43| ) . For larger values of Ay, the 
square root picks the second branch definition. Consequently, s\ and S2 are switched. 
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unnormalized form as 



2iir 



ei 



e 3 



2l7T 

4Ay-2i7r+eT _ 


(S1+S2) 


3J7T 




2l7T 

— 2 A j/+i7r+e a - 


(si+S2) 


1 




2jvr 




4Ay-2i7r+^si- 


4l7T \ 



e 2 



2i?r 



3?7r 



7 4i7T 

-2Ay+i7r+l e~5~ i 
1 



i+«2 



3«7T 



2Ai/+i7r+l s±+e~5~ S2 



V 



(44) 



They depend only on the geometrical parameter Ay. It can be easily checked that in 
the limit Ay — > 00 they become pure color states, t\ a singlet, e 2 and e 3 antisymmetric 
and symmetric combinations of the two color octets in the basis, Eq. (p3|). However, 
for the values of Ay we are interested in, they are color-mixed states. Following Ref. 
|jlj], we will refer to t\ as the quasi- singlet, and to e 2 and e 3 as the two quasi-octets. 
Numerical values of the real and imaginary parts of the eigenvalues are represented in 
Fig. According to the discussion at the end of Sec. EO , the quasi-octet components 
of the hard scattering, for large values of Ay, will have strong suppression factors, 
compared to the quasi- singlet, because their real parts grow much faster. In fact, 
we will see in Sec. |6] that, in this situation, only the quasi-singlet component of the 
scattering survives. 



5.2 Soft anomalous dimension for gg 



99 



We work in the basis of Eq. (p7[) . The structure of the soft anomalous dimension 
matrix is block-diagonal: 



T 5 (At/, Ay) 



T3X3 O3X5 
05x3 T5x5 



(45) 



with the matrix r 3X 3 given by: 



3x3 






-3(6 - P2) 







3p 2 + 327T 



(46) 



and the matrix r 5X 5 given by: 



/ 3^ + p 2 + 2m 




r5x5 — — 



ft, 
2n 










3p 2 
— 3iir 



12 „• 

--=-Z7T 

5 



-12mt 

— 3«7T 

3p 2 


-fiTT 





— 6i7T 



-3(6 - P2) 

— t«7T 








9 • 

5 



3p 2 - 56 / 



(47) 
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The functions £2 and p 2 are very similar to £1 and pi of Eq. 
£ 2 (Ay) = ztt - 2Ay, 

tanh(%2) +tanh(^) 



p 2 (Ar?, Ay) = 2 In 
The eigenvalues of the matrix are 



k tanh(^) - tanh 



+ 3Z7T . 



(48) 
(49) 



2 / 



Ai = 


A 5 


= 3^ 
2vr 


{p2 - m) 


A 2 = 


A 4 


a s 
= 3— 
2tt 


(P2 - 6) 


A 3 = 


A 6 


a s 
= 3— 
2tt 


(P2 + wr) 


A 7 = 


a s 
2tt 


(3p 2 - 


£2 - 4r/ 2 


A 8 = 


a s 
2n~ 


(3p 2 - 


6 + 4^ 2 



(50) 



where we have introduced the following function of the interjet rapidity: 

2 J Ay 2 - in Ay 



V2(Ay) = 2V Ay 2 - mAy - n 2 . (51) 

Due to the block structure of the soft anomalous dimension matrix, the eigenvectors 
corresponding to the first three eigenvalues coincide with ci, c 2 and C3 of Eq. (|27|), 
while the remaining eigenvectors have zero components along these color directions. 
Correspondingly, the structure of the matrix , and consequently of Eqs. 

(|16D and (|i~7D , is block- diagonal, 



-1 







3x5 



(52) 



13x3 

\ 5X 3 ^ '5x5 , 

Referring to the discussion below Eq. fl3~l|) , and to Eqs. (|i"6|) and (|T7|), we see that 
the dimensionality of the hard scattering remains 5x5 after the change from the 
original basis to the basis of the eigenvectors of Tg, because the soft dynamics of the 
physical and unphysical color modes are completely decoupled. Here, once again, the 
eigenvectors depend only on the gap width Ay. Explicitly, we find 

( \ 

• 1 9 o I 0(3) ' 

5 i3 ' 1 ~ ' ' ' Ci ~ I e ^ 5) 

V o( 5 ) J 



4. ..8. 



(53) 



For e^, ef 1 and e 6 5 ^ we have the simple results 

/ -15 
f Of - 2Ay 

4 5) = , ef= f(2m-2Ay) 

3 
1 



3iAy 
2tt 



,(5) 



/ -15 \ 

-3 



(54) 



21 




Figure 7: Plot of the real (left) and imaginary (right) part of the eigenvalues (Aj = 
\/a s ) of the soft anomalous dimension matrix for g + g — > g + g scattering. In both 
plots the solid line identifies the quasi-singlet eigenvalue, A 7 , and the dotted line the 
quasi-27, A 8 . In the plot of the real parts the dashed line identifies the two degenerate 
quasi-octets, A 5 and A 6 , and the short dashed line the quasi 10 + 10. In the plot of 
the imaginary parts the dashed line corresponds to both the quasi-octet A5 and the 
quasi 10 + 10, whereas the dot dashed line identifies the other quasi-octet, Xq. 



It is easy to see that for large interjet regions (Ay — > 00) ef^ is oriented along the 10 x 
10 color component, while and are respectively symmetric and antisymmetric 
combinations of the 85 and 8a color components. The last two eigenvectors have a 
more complicated analytic form: 



with: 



= e 



= e 



(5) 
+ 

(5) 



(55) 



97 

(ei 5) ) i = — — [32Ay 4 - lQiAy 3 (4tt ± irj 2 ) - AnAy 2 (19tt ± 6irj 2 ) 



A/±7T 2 

+4n 2 iAy (IItt ± 5ir] 2 ) + tt 3 (13tt ± 6*772) 

27 r 

9 

X±TT 

-7T 2 (ll7T±7^ 2 ) 



-AAy 2 - 2Ay (-2m ± r] 2 ) + vr (3tt ± ir] 2 
24iAy 3 + 12Ay 2 (3vr ± ir] 2 ) + 2nAy (±6r] : 



17m) 



(5) 



-9tt 

1. 



(2iAy + tt ± 2ir) 2 ) 



(56) 
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In the denominators, the normalization A/± is given by 

M± = -28Ay 2 + AiAy (7tt ± ir] 2 ) + tt (31tt ± 2ir] 2 ) . (57) 

It is easily checked that e\ , in the limit of a very wide interjet region, reduces to a 
pure color singlet state, while eg 5 '' points to the 27 color direction. Again, for typical 
values of Ay, these eigenvectors are not pure color states, but rather reflect mixing 
between the different color components of the hard scattering. From the eigenvalues 
of Eq. (|50D , plotted in Fig. [7], we can see, as already for q + g — > q + g, that the 
quasi-singlet component of the scattering, e 7 , dominates the others in the limit of 
large gaps. 



5.3 Soft anomalous dimension for qq — > qq 

In the basis of Eq. (|3^) , the soft anomalous dimension matrix is 



T s (Ar ] ,Ay) = ^- 



Cp_ 

1 N c 

where the functions £ an d Po are the analogs of those in Eqs. ( |40| ) and ( f£9] 

N?-2 



(58) 



UAy) = -2N e Ay + 2m- 



N 2 - 1 /tanh (AR) + tanh f ^ 
Po(A7 ? ,A 2/ ) = 2^^1n' 



W c \tanh (Aa) - tanh 



2 iv 2 -2 

H Ay - 2«7T — ^ 

Nr. y Nr. 



(59) 



In particular, p$, which depends on the scattering angle, differs slightly from the 
corresponding function computed in Ref. [IJ for the Tevatron gaps, because of the 
different underlying kinematics. On the other hand £o> depending only on the geom- 
etry, through Ay, is identical. The two eigenvalues are 



Ai (At;, Ay) 



2tt 



-On == 

2 H 2JNZ 



N c e - 32C F 7T 2 



X 2 (A V ,Ay) = -i 
corresponding to the eigenvectors 



2 P0+ 2VWr 



1 



N c $ - 32C f tt< 



(60) 



ei 



('2 



I 2ZL 

V 1 



x 6 



/iVc 



N c {UAy))-32C F ir< 



±[Co + ^N c (UAy)Y-32C F ^ 
1 



(61) 
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Figure 8: Plot of the real (left) and imaginary (right) part of the eigenvalues (Aj = 
Xi/a s ) of the soft anomalous dimension matrix for q + q —>■ q + q scattering. The solid 
line identifies the quasi-singlet eigenvalue, Ai, the short dashed line the quasi-octet, 
A 2 . 



which coincide precisely with the ones given in Ref. Jl|, since they depend only on 
the geometric parameter Ay. In the limit of an asymptotically large gap, e\ becomes 
a pure color singlet vector, while e 2 orients itself along the octet direction. Referring 
to ei and e 2 as quasi-singlet and quasi-octet, and looking at the behavior of the real 
parts of the eigenvalues, Fig. [5], we can predict that, also for this partonic process, the 
quasi-singlet component of the scattering in Eq. ( |14"D will dominate the quasi-octet in 
the large gap limit. 



5.4 Soft anomalous dimension for qq —> qq 

In the basis of Eq. (|3^) , the soft anomalous dimension matrix is 



T 5 ( A V , Ay) 



Air 



PO + & ~ W c 4f Z7T 



Ai-K I ) 
N c J 



where the functions p' and £q are: 



e (Ay) = -2N c Ay + 



p' (A V ,Ay) = 2 



tanh ($p) + tanh f ^ 

Id I f^i hf4 I 2/tt 

tanh (^) - tanh 



H Ay 



(62) 



(63) 



The two eigenvalues are: 
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Re (A) Im(A) 




Figure 9: Plot of the real (left) and imaginary (right) part of the eigenvalues (Aj = 
Xi/a s ) of the soft anomalous dimension matrix for q + q —>■ q + q scattering. The solid 
line identifies the quasi-singlet eigenvalue, Ai, the short dashed line the quasi-octet, 
A 2 . 



Ai (A77, Ay) 



2^ 



A2 (A77, Ay) = g 



1 , 1 

2 Po 2JW r 



k+ 1 — 



iV^ 2 - 32^C F - ^ 
2i?r' 



— v - 32tt 2 C f - — 



which correspond to the eigenvectors: 



ei 



e 2 



47rCf / 



V 



- v & - 32tt 2 C f 
1 



AttCf 




(64) 



(65) 



Notice the similarity with the process qq — > qq. We can again identify in e\ the quasi- 
singlet, and in e 2 the quasi-octet. For large gaps, the quasi-octet component of the 
scattering will be strongly suppressed with respect to the quasi-singlet, because the 
real part of the corresponding eigenvalue grows much faster with Arj, as graphically 
shown in Fig. ||. 



5.5 Soft anomalous dimension for 7$ 



qq 



For the direct process 75 — > qq, which, from Fig. ||, gives the second largest con- 
tribution to the denominator of the fraction, after the partonic reaction qg — > qg, 
there is only one possible color flow at the hard scattering. Thus, the soft anomalous 
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Figure 10: Plot of the real (left) and imaginary (right) part of the soft anomalous 
dimension for 7 + g — > q + q scattering. 



dimension reduces to a function, given by 

'tanh (fz) + tanh 



ft* 



T s (A V , Ay) = ^ 



C F ln 



tanh (^2) - tanh (^) 



+ 2mC A 



(66) 



whose real and imaginary [] parts are plotted in Fig. [It]. We see that the real part 
of Ts grows faster with At] than the real parts of the quasi-singlet eigenvalues for all 
the resolved processes, apart from gg — > gg, Figs. ^ - In fact, we will see in the 
next section that the corresponding contribution to the gap cross section will have a 
fast decay rate. 



6 Numerical Results and Discussion 

From the previous two sections, we have all the tools necessary for the evaluation of the 
partonic gap cross section, Eq. (|Tj|), which we convolute with the parton distributions 
and with the photon distribution in the electron, according to Eq. ([7|), to get the full 
gap cross section. The necessary numerical integrals have been performed with the 
routine VEGAS. We have used the CTEQ4L parton distribution set |2?J for the 



proton, and the GRV-G LO set |27J for the photon. For the interjet energy threshold 
Qo, which in our approach defines the gap event (see Sec. ||), we have considered 
several different values, from Q = 350MeV up to Qq = 1050MeV, and we have 
evaluated the related gap cross section as a function of the pseudorapidity difference 
between the jets, At]. Figs. |TT|a,b show our results for the gap cross section and 
the corresponding gap fraction. Both quantities increase with increasing Q , simply 



3 The imaginary part of Eq. ( |66| ) can be extracted directly from Ref. [p5| . 
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because a stronger limitation on interjet radiation leads to a stronger suppression 
of the result. For the most forward configuration of the dijet analyzed in Ref. ||, 



A 77 = 3.75, the gap fraction, Fig. |TT]b, varies from about 10% at Q = 350MeV to 
about 25% at Q = 1050MeV. One may ask, for both the gap cross section and the 
gap fraction, whether any of the plots in Figs. pJ]a,b fits the experimental data of 
Ref. 0, for which, as discussed in Sec. [| the gap event is identified in a related, but 
different way. With the choice Qq = 350MeV for the interjet energy, we obtain a good 



match, as shown in Fig. Oa,b. Observe that, in this situation, Q is not much larger 



than A in Eq. ([18]). Results in this region, at the boundary with non-perturbative 
QCD, can be considered as smooth extrapolations of perturbative QCD resummation. 



In Fig. [13], the analog of Fig. 3d of Ref. ||, we show again our gap fraction at 
Qo = 350MeV, from Eqs. (|1^) and (|T9"D, and compare it with the result of a two 
parameter fit to the expression 

f\a 1 ^Ar 1 ) = C{a 1 P)e a ^ + p. (67) 
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as done in Ref. ||. We find a = —2.1 and (3 = 8.3%, in good agreement with the 
values of Ref. 0: a = —2.7 ± 0.3 and (3 = (7 ± 2)%. Here, in analogy with Ref. 
||, we have imposed the constraint / fit = 1 at Arj = 2, meaning that the gap cross 
section has to reduce to the full dijet cross section in the absence of a gap. Notice 
that this information is not encoded in Eq. ([Tj|), because the last two factors on the 
right hand side do not reduce exactly to unity for Arj = 2. This can be seen from 
Eq. ( |15|) and from the eigenvalues of the soft anomalous dimension matrices, Figs. 
|6|-|1~0"|. In the limit Arj — > 2, the real parts of the eigenvalues vanish, but the imaginary 
parts do not, causing a slight suppression in the mixed terms, {(3 ^ 7), of the hard 
scattering matrix in Eq. ([18]). The validity of our formula for the gap cross section, 



Eq. (|18D, is restricted to interjet rapidity widths Ay neither too small, nor too large. 
We will limit ourselves to the Arj range 2.6 < At] < 4.0. 

The reasons for these restrictions are as follows. Our expression in Eq. (|i~8| ) is the 

result of the resummation of logarithms of the soft interjet energy, a" (— i) ln n (<5o/ \f—i)-, 
which are assumed to be the only large, perturbatively dangerous quantities in the 
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Figure 13: Analog of Fig. 3d of Ref. 0. The gap fraction of Fig. |T2| (solid line) is 
redisplayed and compared with the result of a fit to an exponential plus a constant. 



problem. When Ay — > 0, the gap closes and the cross section becomes inclusive. 

More specifically, logarithms of would be suppressed by powers of Ay. Analo- 
gously, when Ay — > oo (i. e. , At] — > oo), terms of the form a£(— t)(Ay) n , coming 
from the eigenvalues of the Fg s, are large at each order in perturbation theory, and 
we need further resummations. We should also emphasize that the limit At] — > oo 
corresponds to the Regge region of the scattering, when s becomes large at fixed 
momentum transfer, t. Resummation in the Regge limit organizes logarithms of the 
form a"ln n (s/(— coming from the BFKL ladders of gluons ||. For double loga- 
rithmic terms like a™ ln n (s/ (—ij) \n n (Q / y—i), which occur in Reggeized color octet 
exchange, both BFKL resummation and our method give the same result. 

Note that the suppression we find here is not double logarithmic in Qq |28] . Rather, 
it comes from the exponentiation of single "soft" logarithms only, as shown in Eq. 
(|IB|). The underlying physical reason is that a soft gluon emitted into the interjet 
region can never become collinear to either of the forward partons from which the 
jets originate. 

Finally, in Figs. p^ - |T7| we show the contributions to the gap cross section, Eq. 
(|18|) at Qq = 350MeV, from the different partonic processes studied in Sec. |[ In 
each case we give the full result and its decomposition into quasi-color components. 
We see that, as anticipated in Sec. [|, the quasi-singlet components, in all cases, give 
almost the whole gap cross section at large values of Ar/. It then becomes clear that 
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Figure 14: Overall contribution to the gap cross section from qg — > qg scattering 
(solid line). Contribution from the quasi-singlet component (dashed) and from the 
two quasi-octet components (dotted and dot dashed). 



the parameter (3 = 8.3%, extracted from the fit, Eq. fl67|), can be associated with 
the asymptotic quasi-singlet fraction. From Figs. ^|-^ we see that at high Ar] the real 
parts of the quasi-singlet eigenvalues are ordered, from the lowest value of q q — > q q 
scattering, through the intermediate values of processes qg — > qg and q q — > q q, up to 
the largest value of g g — ► g g scattering. In Figs. pT4| ~ pT7| , the decay slopes of the quasi- 
singlet components of the corresponding contributions to the gap cross section reflect 
the same ordering property, as a consequence of Eq. (|I4"D . This observation indicates 
that gaps tend to be more easily formed with scattered quarks and antiquarks than 
with gluons |J. 



7 Conclusions 

In this paper, we have shown that it is possible to analyze the dijet rapidity gap events 
from photoproduction, observed at HERA, by introducing an energy- dependent def- 
inition of the gap, as already done in Ref. |lj for the Tevatron events. The result we 
have found is perturbative, as the ordering of the different scales in Eq. ([14]) shows, 

with A < Q c < y—i. The experimental behavior of the gap fraction is approxi- 
mately reproduced by fixing the threshold of interjet energy, which defines the gap, 
at Qq = 350MeV. Our result also predicts how gap fractions increase with Q . We 
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Figure 15: Overall contribution to the gap cross section from gg — > gg scattering 
(solid line). Contribution from the quasi-singlet component (dashed) and from the 
two quasi-octet components (dotted and dot dashed). The contributions of the quasi- 
10 x 10 and of the quasi-27 component are not exhibited, because they are too small. 



conjecture that in photoproduction non-perturbative "survival" effects from the in- 
teractions of spectator partons, which can give rise to interjet multiplicity |T7[], should 
be reduced with respect to the case of pp scattering, because there is only one incom- 
ing hadron. The experimental determination of the double differential cross section 
d 2 a ep /dAr]dQ c , (and related fraction), Eq. (|18"D, for different values of the interjet 
energy flow identifying the gap, as in Fig. |ll], would offer a significant test of the 
perturbative dynamics of QCD radiation. 
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Figure 16: Overall contribution to the gap cross section from qq — > qq scattering 
(solid line). Contribution from the quasi-singlet component (dashed) and from the 
quasi-octet component (dotted). 



A The Denominator of the Gap Fraction: Par- 
tonic Cross Sections and Parton Luminosities 

In this Appendix we will summarize for completeness the Born cross sections, and 
corresponding parton luminosities, for all the partonic processes contributing 
to the denominator of the gap fraction, Eq. (|19|). We will consider for the quarks 
four active flavors, nj = 4, and assume flavor symmetry. Correspondingly, the quark 
luminosities will be expressed as functions of the proton valence distributions, 4> u / p 
and 4>d/ P , and of the sea distributions for both the proton and the photon, respectively 
(fisea/p and sea / 7 . The analytic cross sections will always be of the form 



dAr] ~ V V 2 J J ~ sw V2s 



2c osh 2 (^L)) 1 a 2 s (t) (£) \M^(s,t,u)\ 2 



2 



'Ar] 




where, in both equalities, the first factor is the Jacobian [dcos(9)/dArj], and the rest 
is the partonic cross section da^ / d cos(#) . The second equality defines the normal- 
ization of the hard scattering matrix in Eq. (|TT|) . in the list below, we will specify 
the invariant matrix elements squared, \A4^\s, i, u)\ 2 [ I4|, p4 . 
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Figure 17: Overall contribution to the gap cross section from qq — > qq scattering 
(solid line). Contribution from the quasi-singlet component (dashed) and from the 
quasi-octet component (dotted). Negative interference terms [l], corresponding to 
(3 7^ 7 in Eq. flUf ), are not exhibited separately. 



Processes g(g)# -> g(g)#, flrg(g) 



^ = $9h (fiu/p + + 60 sea /p) + 80 sea / 7 9 /p 



+ s 2 4 s 2 + u 2 



l-M (f) (s,t,^)| 2 = q ~ • 69 

t 2 9 SM 

Processes q(q)g — > gq(q), gq(q) — > q{<i)g (observe that they differ from the 
processes of Eq. (|S9|) by the exchange of particles 1 and 2 in Eq. fl2"T|) , and, 
correspondingly, of the Mandelstam invariants t and u)\ 



sea/p 

t 2 + s 2 4 s 2 + t 2 
m 2 9 si 



\MV(s,t,u)\ 2 = L+L _££__+! (70) 



Process — > 

L (f) = 9 / 7 9 / p 



. ,m,. - -m9 9 / su tu st\ ,„ . 

I^ (£) (*.*.«)l a = 2 3 "F + i^ + ^ ■ (71) 
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Process gg — ► qq: 



[(f>u/p + 4>d/p + 60sca/p) 



\M^(s,t,u)\ 2 = - 



Process qq — > gg : 



|-M (f) (U«)| 2 
Process gg' — > q'q: 

L® 

Processes qq — > qq, qq — > gg: 

L<*> = 

\M^(s,t,u)\ 2 = 

Processes gg' — > gg', gg' — > gg': 

L<*> 

i-M (f) (M,^r 

Processes gg — > g'g', gg — > g'g': 

L<*> 

|A* (f) (s,f,«)| 2 

Processes gg — ► gg, gg — ► gg: 

L« = 

|7W (f) (S,t,M)| 2 = 

Processes gg' — > g'g, gg' — > g'g: 

L<*> 

|A* (f) (S,t,fi)| 2 



4 s 2 + u 2 A s 2 + t 2 8 s 2 
9 £2 +g ^2 ^7£^ 



'sea/p 



3</>sea/7 (0u/p + 

4 S 2 + u 2 

9 



30sea/7 (^u/p + 0d/p + 60 sca / p ) 
4 5 2 +t 2 
9 W 2 ' 



0sea/ 7 {4>u/p + 4>d/p + 60 

sca/p 

4 s 2 + u 2 4 t 2 + m 2 8 m 2 
9 ~T 2 + 9 ^ 27 it 



30 S ea/ 7 [4>u/p + <Pd/p + 60 sea/p 

4 s 2 + M 2 

9 £ 



30sea/7 {^Pu/p 

At 2 + u 2 



9 s 2 



0sca/7 (^0u/p + 0d/p + 60 

sca/p 

4 s 2 + g 4 P + n 2 8 g 
9 w 2 + 9 s 2 27 



30sea/7 (0u/p + 0d/p + 60 sca / p ) 

4 g 2 +t 2 
9 M 2 
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Processes gg — > q'q', qq — » q'q': 

L (f) = 30 sca/7 (<f) u/p + (f) d/p + 60 sca/p ) 

\M^\s,i,u)\ 2 = (80) 
Processes qq — ^ gg, qq — > 

(lU) |. = --_ + -_. (8!) 
Processes gg — > gg, (7(7 — * gg 

L (f) = <t> gh <t> g /p 



Processes 7g — > gg, 75 — > gg 



10 a ( 

9 a s \i u 

Process 7g(g) — > g(g)g, 



M^(s,t,u)\ 2 = ^-^(t + t). (83) 



r(H 4 _, ^ 15 x 

^ = gPu/p + gPd/p + — 0sca/p 

l-M«(5,t,n)| 2 = ^ f— + — ) • (84) 
a s \ s -u) 
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